Dynamics of Localized Waves 
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We have measured pulsed microwave transmission through quasi-lD samples with lengths up to 
three localization lengths. For times approaching four times the diffusion time td, transmission 
is diffusive in accord with the self-consistent theory of localization for the renormalized diffusion 
coefficient in space and frequency, D(z, St). For longer times, the transmission decay rate first agrees 
with and later falls increasingly below the self-consistent theory. Beyond the Heisenberg time, the 
decay rate approaches the predictions of a dynamic single parameter scaling model which reflects 
the decay of long-lived localized modes and converges to the results of ID simulations. 
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The theory of localization was developed in the context 
of electronic conduction and has been widely applied to 
steady state transport [l|, |2| . Because of electron-electron 
interactions, however, the localization transition is not a 
pure single-particle Anderson transition. The descrip- 
tion of localization is made all the more complex by in- 
elastic scattering, since the impact of scattering grows 
with pathlength and waves following paths of all lengths 
contribute to the conductance. Localization can also be 
studied for classical waves in optics and acoustics 
However, the exponential decay of transmission found 
in the presence of absorption even for diffusive waves 
makes it difficult to determine the localization length 
from measurements of the exponential scaling of trans- 
mission [J, |5|, [6[ . Similarly, the localization length cannot 
be determined directly from the rounding of the coherent 
backscattering peak which is produced by both absorp- 
tion and localization [H, 0| . 

The ensemble average of pulsed transmission, (I(t)), 
is of particular interest for classical waves [3] because it 
allows the increasing impact of weak localization to be 
disentangled from absorption, since the relative weights 
of paths within the sample at a g iven delay time is unaf- 
fected by absorption [1, 0, [13, [Tl( . Far from the localiza- 
tion threshold, measurements of (I(t)} in opaque samples 
have generally been well described by diffusion theory 0] . 
For times greater than the diffusion time, rrj, higher dif- 
fusion modes decay rapidly leaving energy in the lowest 
diffusion mode so that the decay rate approaches the con- 
stant, 1/t d = 7r 2 D I (L+2z ) 2 0. Here, D is the diffusion 
coefficient, L is the sample length, and Zq is the distance 
beyond the boundary at which the intensity within the 
sample extrapolates to zero 12j. However, a progressive 
suppression of the decay rate has been observed in recent 
microwave [icl | and optical fl~0 ] measurements in strongly 
scattering samples in which steady state transmission is 
essentially diffusive. 

The impact of localization on electron dynamics for 
single electrons at T=Q had been calculated using dia- 



grammatic, nonlinear a, and supersymmetry approaches 
[13j |. Localization is achieved when the average spac- 
ing between quasimodes exceeds their average linewidth, 
Au>5v 1A\. In the time domain, this is the condi- 



tion that the Thouless time exceeds the Heisenberg time, 
TTh>TH, where TTi i =l/5v=ir 2 Tr) 1 and th=1/Ai/ is the 
time required to visit each coherence volume of the sam- 
ple. 

The slowing decay of (I{t)) reflects the increasing 
prominence of longer-lived modes which are more re- 
mote from the sample boundaries or are more sharply 
peaked within the sample 

S E 0, E3- This is as- 
sociated with the increasing enhancement of weak lo- 
calization with longer pathlength due to the scattering 
of the wave as it crosses over its trajectory. Vollhardt 
and Wolfle (VW) developed a self-consistent diagram- 
matic theory of localization within a medium in terms of 
a frequency-dependent renormalized diffusion coefficient, 
D{Q) [18]. To be self-consistent, Van Tiggelen et al. [ljS| 
argued that D must also be a function of depth within 
a bounded sample, D(z, f2). Skipetrov and Van Tiggelen 
used the self-consistent localization theory (SCLT) to de- 
scribe waves near the mobility edge for t<Tn in quas i- ID 
(20| and slab geometries [21|. They described [2fj | key 
features observed in microwave measurements for diffu- 
sive waves in quasi-lD [l(| and found a l/t 2 falloff in 
reflection for localized waves. 

In this Letter, we present microwave measurements 
of dynamic transmission for localized waves in quasi- 
lD samples with L greater than the average localiza- 
tion length, £. Four different approaches have been 
used to analyze the measurements in different time 
ranges. For times up to several times the peak ar- 
rival times, tp^Tr), (I(t)) can be well described by 
a simple diffusion theory. These results are in ac- 
cord with the SCLT, which includes a position and fre- 
quency dependent diffusion coefficient D(z,Q) [20L 21 1, 



and suggest that the renormalization of D is insignifi- 
cant at early times. For £>4td, the transmission decay 
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rate, A(t)=-(d(I(t))/dt)/(I(t))=-d\n(I(t))/dt, is pro- 
gressively suppressed by localization up to a factor of 
nearly 3 relative to the early diffusive decay rate, 1/td. 
Measurements of A(t) are compared to the SCLT, to 
ID simulations, and to a dynamic single parameter scal- 
ing (SPS) model. This model is based upon a Gaus- 
sian distribution of Lyapunov exponents, 7=l/2£;, with 
var(7)=7/L [23]. Self-consistent calculations provide 
reasonable agreement with measurements for KAtd, but 
give substantially higher values for A(t) at longer times. 
ID simulations give a peak in the decay rate which is 
higher and peaks later than measurements. The dy- 
namic SPS model rises within the pulse width and falls 
below measurements for t<rn. The diffusion-like delay of 
the transmission peak reflects the impact of short-lived, 
spectrally-overlapping, quasi-extende d q uasimodes, de- 
scribed by Pendry as necklace states [l5j. These results 
indicate the greater prominence of overlapping modes in 
ID than in quasi-lD samples. For t>Ta, the results of 
ID simulations and the dynamic SPS model which reflect 
the contributions of long-lived modes with £<£ converge 
and are in good agreement with measurements. 

Microwave spectra of the field transmitted through 
low-density random mixtures of alumina spheres were 
taken with the use of a vector network analyzer. The 
wave is launched and detected by conical horns placed 30 
cm in front of and behind the sample. Alumina spheres 
with diameter 0.95 cm and index of refraction 3.14 are 
embedded within Styrofoam shells to produce a sample 
with alumina volume fraction 0.068 which displays dis- 
tinct sphere resonances [23j |. The sample is contained 
within a copper tube with diameter of 7.3 cm and plas- 
tic end pieces. Spectra are taken for 10,000 configura- 
tions produced by briefly rotating the tube, in samples 
of length 13, 29, 40, 50, 61, and 90 cm. Measurements 
are made just above the first sphere resonance over the 
frequency range 9.95-10.15 GHz, in which the change in 
static and dynamic propagation parameters is small. 

The time response to a Gaussian intensity pulse of 
width (2\/27rcr) _1 peaked at t=0 was obtained by taking 
the Fourier transform of the field spectrum multiplied by 
a Gaussian envelope of width a centered at f c . The field 
of the time response is squared to give the transmitted 
intensity I(t) for each sample realization. The average 
transmitted intensity (I(t)} is found by averaging over 
the ensemble, and then over the frequency interval by 
shifting f c . The measured time response includes a con- 
stant background at long times [10(, which results from 
noise in the field spectra. Subtracting this background 
enhanced the dynamic range by 16 dB. 

To compensate for losses due to absorption, and 
thus to facilitate the comparison of the measurements 
to dynamical models of localization, (/(<)) is multi- 
plied by exp(i/r a ), where l/r a is the absorption rate. 
1/t o =Q.0064 ns" 1 is found from the decay rate of trans- 
mission in a 40cm-long sample with copper end caps, 
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FIG. 1: (color on line) (a) The transmission decay rate in a 
L=40 cm sample weakly coupled to the measurement ports, 
so that the leakage rate is well below the absorption rate l/r . 
l/r o =0. 0064 ns -1 is found as the average of A(t) for t>1.2 fis, 
multiplied by 0.91 to correct for the copper end caps added 
to the tube; (b) The average time response to a Gaussian 
pulse of cr=15 MHz in the L=61 cm sample is compared to 
the SCLT and diffusion theory. All the curves are normalized 
to be unity at the peak. The inset shows a fit of the SCLT to 
the measured data at early times. 



which is weakly coupled to the measurement ports, so 
that the leakage rate is well below the absorption rate 
[Fig. la]. 

Measurements are compared to calculations of the 
renormalized diffusion constant, D(z, Q). In an open sys- 
tem, the SCLT of VW [l8j] can be generalized as follows 
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where Db=ve(-/3 is the Boltzmann diffusion constant, 
I is the transport mean free path, 7V=32 is the number 
of transverse propagating channels in the sample, and 
Ue=11.85 cm/ns is the transport velocity at 10 GHz [25| . 
The diagonal intensity Green function, G{z, z; repre- 
sents the return probability at z and can be obtained 
from the following generalized diffusion equation, 

d z [D(z,n)d z G(z,z';Cl)] + iClG(z,z';Q) = -S(z-z'), 

(2) 

with mixed boundary conditions at both ends, 
ZoD(zi,; Vt)d z G{zb, z'\ fi) T D B G(zb, z'\ £1)=0, where z b =Q 
or L, zo=(2/3)£(l+R)/(l—R) is the extrapolation 
length, and R is the internal reflection coefficient. Equa- 
tions (1) and (2) are solved self-consistently in real f2- 
space to obtain G(z, z'; f2). The intensity just outside the 
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FIG. 2: (color on line) Measured transmission peak arrival 
time in alumina samples (solid circles), following incident 
Gaussian pulse with ct-Ct^ 1 , and the prediction of diffusion 
theory, t p = 0.89ro, with 73=7.9 cm 2 /ns and 2o=6.1 cm (solid 
line) . 



output surface, (I(t)), is obtained by taking the Fourier 
transform of G(L, z'=£; f2) in fl. A fit of the expression 
for (I(t)) at early times to the measured data for sam- 
ples with L=6l cm, with £ and R as fitting parameters, 
gives £=2.0 cm and i?=0.64, and is shown in the inset 
of Fig. 1(b). Excellent agreements between the SCLT 
and experiment is found for t<4r£>. To compare the 
long-time behavior, we present these results in a semi- 
log plot in Fig. 1(b). Also shown in Fig. 1(b) is the 
result of diffusion theory in which localization effects are 
absent. The surprising agreement between diffusion the- 
ory, SCLT, and measurement for t<2ro suggests that 
the renormalization of the diffusion constant is negligi- 
ble at early times even for localized waves studied here. 
The applicability of diffusion theory at early times in our 
samples can also be seen in the excellent agreement be- 
tween the measured and calculated peak arrival times, t p , 
shown in Fig. 2. We believe, this is because TD<TH<7Th 
in our samples. The first inequality guarantees that the 
wave propagation is diffusive for early times, while the 
second indicates that the samples are in the localized 
regime. For example, we have td=70.4 ns, th=528 ns, 
and Txh=695 ns for L=61cm and td=136.3 ns, th=780 
ns, and TTh=1346 ns for _L=90cm, respectively. The val- 
ues of th are obtained from measurements of the average 
spacing between modes in a closed sample [23j]. 

Since the SCLT includes localization effects, it gives a 
slower-than-exponential decay in (I(t)} seen in Fig. 1(b). 
However, for i>4T£>, the SCLT underestimates the local- 
ization effects and thus gives a higher decay rate than 
measured. To show this more clearly, we plot both the 
measured and calculated decay rates in Fig. 3. It is also 
seen that, at a given delay time, a larger deviation is 
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FIG. 3: (color on line) The transmission decay rate, A(t), 
measured in samples with L=61 (a) and 90 cm (b) is com- 
pared to predictions of the SCLT and dynamic SPS model. 
The comparison to ID simulations is made by using the di- 
mensionless decay rate, tA, plotted as a function of the di- 
mensionless time, t/r, where t=£/ /3ve=50.6 ns. 



found for the L=61 cm sample than for the L=90 cm 
sample which is deeper in the localization regime. This 
is because, for a given delay, the ratio of pathlength to 
sample volume is greater in the shorter sample, so that 
the number of closed loops and the consequent renormal- 
ization is greater in the shorter sample. A slower decay 
in both samples suggests that more extended quasimodcs 
that may exhibit diffusion-like behavior, and which can 
be described by the SCLT, have largely decayed so that 
the wave is transmitted increasingly via long-lived local- 
ized modes. 

Because transmission at long times may be determined 
by the longitudinal structure of spectrally isolated local- 
ized modes, we consider the statistics of such modes. In 
ID, the steady-state intensity at z=L of resonantly ex- 
cited modes relative to the incident wave at 2=0+, for 
modes peaked a distance z from either sample boundary, 
was given by Azbel [3, T=exp(-27(L - 2z)). The Lya- 
punov exponent, 7=l/2£, is drawn from a Gaussian dis- 
tribution, P("/) = \/L/2itj exp[— (7— 7) 2 /(27/L)1, accord- 
ing to the SPS hypothesis [22[. We assume the position 
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of peak intensity for the modes is uniformly distributed 
between and L. The decay rate of localized states is 
the ratio of the sum of the outgoing flux at the open ends 
to the integrated wave energy inside the sample, 



r(7, z ) = Pv E -, 



1 + exp(-27(L - 2z)) 



[2 exp(2z 7 ) - exp(-27(L - 2z)) - l]/2 7 ' 

(3) 

The coupling factor, 0, reflects the reduction of the trans- 
mitted flux due to the angular spread about the normal 
to the interface as well as the angular average of inter- 
nal reflection at the interface, R, and by the character of 
transport near the boundary. The transmitted intensity 
is then, 



(/(*)> - 



1 

2L 



e?7 
2/L Jo 



L/2 



dzP(j)Tr z exp(~Tt) . (4) 



The factor T (T/2) 2 exp(— Ft) is the square of the Fourier 
transform of a Lorentzian line in the field spectrum with 
linewidth T/2ir. Internal reflectivity will not change T 
appreciably since the intensity throughout the sample is 
enhanced by the same factor that inhibits transmission 
at the interface. By using £=30 cm (23|, we fit A(t) calcu- 
lated from Eq. (4) to measurements in the samples with 
L=61 and 90 cm, to find /3=0.05. The small coupling fac- 
tor at the boundary is a consequence both of the strong 
reflection and the suppressed flow of energy in the expo- 
nential tail of localized modes. These results are shown in 
Fig. 3. The rise time of transmission in this model is the 
rise time of a localized mode, which is essentially the rise 
time of the incident pulse. Thus the slower rise observed 
in the experiment indicates the dominance of overlapping 
as opposed to isolated modes. At long times, £>th, how- 
ever, the excellent agreement of the dynamic SPS model 
with measurements indicates that the energy within the 
sample is stored in long-lived localized modes. 

The above calculations of localized mode dynamics 
may be compared to ID simulations in samples with 
the same values of L/£. We consider a random sample 
of L/a layers with equal thickness, a, embedded in air. 
The dielectric constant in each layer is a random num- 
ber uniformly distributed about e=l, from 0.3 to 1.7. 
A Gaussian pulse with carrier frequency of tOo=1.65c/a, 
where c is the speed of light, and width a— 0.14c/ a is 
incident upon the sample, and the intensity I(t) just be- 
yond the output surface is calculated. Over the width 
of the incident spectrum, £=22a. (1(f)) is obtained by 
averaging over 10,000 configurations. The comparison 
of propagation in quasi-lD and ID samples is facilitated 
by employing a dimcnsionless time t/r, where t—^/(3ve 
[Fig. 3]. The dimensionless decay rate, rA(t), obtained 
from Eq. (4) is a universal function of t/r, which de- 
pends only on the ratio L/£. In ID systems, we take 
(3=1, since the average index of the sample is the same 
as its surroundings and no angle average is necessary. 



Further, we assume ve=c, the effective medium veloc- 
ity, since there are no internal resonances over the pulse 
bandwidth. The greater delay of the peak in tA and 
the higher decay rate after the peak found in ID sim- 
ulations relative to the quasi-lD measurements, as seen 
in Fig. 3, indicate that overlapping quasi-extended modes 
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produce diffusion-like behavior at early 
times are more prominent in ID than in quasi-lD. 

In conclusion, we find that the nature of pulsed trans- 
mission in localized samples evolves with time. At short 
times, t<2ro, propagation is diffusive; at intermediate 
times, 2r£)<t<4ru, transport can be described in terms 
of a position and frequency dependent renormalized dif- 
fusion coefficient, while at later times, £>th, energy flows 
exclusively from isolated localized modes. This transi- 
tion mirrors a change in the distribution of modes over 
time with the distribution flowing from short-lived over- 
lapping modes towards long-lived localized modes with 
increasing time. This work shows that a unified theoret- 
ical description of pulsed propagation of localized waves 
will need to incorporate the full distribution of spacings 
and widths of quasimodes of the random medium as a 
function of the average overlap parameter, 5=Sv/ Av. 
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